Are Galaxy Clusters Suggesting an Accelerating Universe? 
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Abstract. The present cosmic accelerating stage is discussed through a new kinematic method based on the Sunyaev- 
Zel'dovich effect (SZE) and X-ray surface brightness data from galaxy clusters. By using the SZE/X-ray data from 38 galaxy 
clusters in the redshift range 0. 14 < z < 0.89 [Bonamente et al., Astrop. J. 647, 25 (2006)] it is found that the present Universe 
is accelerating and that the transition from an earlier decelerating to a late time accelerating regime is relatively recent. The 
ability of the ongoing Planck satellite mission to obtain tighter constraints on the expansion history through SZE/X-ray angular 
diameters is also discussed. Our results are fully independent on the validity of any metric gravity theory, the possible matter- 
energy contents filling the Universe, as well as on the SNe la Hubble diagram from which the presenting accelerating stage 
was inferred. 
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INTRODUCTION 

One decade ago, observations from distant type la supernovae lead to a landmark conclusion: the universal expansion 
is speeding up and not slowing down as believed since the early days of observational cosmology [1, 2]. This 
phenomenon is normally interpreted as a dynamic influence of some sort of dark energy whose main effect is to 
change the sign of the universal decelerating parameter q(z)[3]. Another possibility is that the cosmic acceleration is a 
manifestation of new gravitational physics (rather than dark energy) that involves a modification of the left hand side 
(geometric sector) of the Einstein field equations. In this kind of theory the Friedmann equation is modified and a late 
time accelerating stage of the Universe is obtained even for a Universe filled only with cold dark matter (CDM). At 
present, the space parameter associated with the cosmic expansion is too degenerate, and as such, it is not possible to 
decode which mechanism or dark energy component is operating in the cosmic dynamics [3, 4]. 

SNe type la are not only the powerful standard candles available but still provides a unique direct access to the 
late time accelerating stage of the Universe. Naturally, this a rather uncomfortable situation from the observational and 
theoretical viewpoints even considering that ten years later, the main observational concerns about errors in SNe type la 
measurements, like host galaxy extinction, intrinsic evolution, possible selection bias in the low redshift sample seem 
to be under control [5]. A promising estimator fully independent of SNe type la and other calibrators of the cosmic 
distance ladder is the angular diameter distance (DA(z)) from a given set of distant objects. It has also been recognized 
that the combination of SZE and X-ray surface brightness measurements may provide useful angular diameters from 
galaxy clusters [6, 7, 8, 9, 10, 1 1]. 

On the other hand, since the mechanism causing the acceleration is still unknown, it is interesting to investigate 
the potentialities of SZE/X-ray technique from a more general viewpoint, that is, independent of the gravity theory 
and the matter-energy contents filling the Universe. The better strategy available so far is to consider the same kind 
of kinematic approach which has been successfully applied for determining the transition deceleration/acceleration in 
the past by using SNe type la measurements [12, 13, 14, 15, 16]. 

In this work, we employ a purely kinematic description of the universal expansion based on angular diameter 
distances of clusters for two different expansions of the deceleration parameter. As we shall see, by using the 
Bonamente et al. [10] sample we find that a kinematic analysis based uniquely on cluster data suggests that the 
Universe undergone a dynamic phase transition (deceleration/acceleration) in a redshift z ~ 0.3. Further, it is also 
shown that the Planck satellite mission data must provide very restrictive limits on the space parameter, thereby 
opening an alternative route for accessing the expansion history of the Universe. 



ANGULAR DIAMETER AND KINEMATIC APPROACH 



Let us now assume that the Universe is spatially flat as motivated by inflation and WMAP measurements [17]. In this 
case, the angular diameter distance in the FRW metric is defined by (in our units c = 1), 



n n+^tf- 1 f du {l+z) ~ l 



rz ru 

/ exp — / [\ + q(u)]d\n(\ + u) 
Jo Jo 



du, 



(1) 



where H(z) = a/a is the Hubble parameter, and, q(z), the deceleration parameter, is defined by 
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In the framework of a flat FRW metric, Eq. (1) is an exact expression for the angular diameter distance. As one 
may check, in the case of a linear two-parameter expansion for q(z) = qo+zqi, the above integral can analytically be 
represented as 
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where qo and q\ are the values of q(z) and its redshift derivative, dq/dz evaluated at z — while 7 is an incomplete 
gamma [18]. By using the above expressions we may get information about qo, q\ and, therefore, about the global 
behavior of q(z). In principle, a dynamic "phase transition" (from decelerating to accelerating) happens at q(z t ) = 0, or 
equivalently, z t = —qo/li- Another interesting parametrization is q{z) = q () + q\z/{l+z) [15, 16]. It has the advantage 
to be well behaved at high redshift while the linear approach diverges at the distant past. Now, the integral (1) assumes 
the form: 

Da(z) = ^^^9r (90+?l) [7(9i +90,9i) 
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where q\ now is the parameter yielding the total correction in the distant past (z ^> 0,q(z) = qo + 1i) an d 7 is again 
the incomplete gamma function. Note that in this case the transition redshift is defined by z t = —qo/{qo + q\). 



SZE/X-RAY TECHNIQUE AND CONSTRAINTS 

The SZE is a small distortion on the Cosmic Microwave Background (CMB) spectrum provoked by the inverse 
Compton scattering of the CMB photons passing through a population of hot electrons. This distortion is characterized 
by a low frequency (< 218 GHz) decrement and higher frequency (> 218 GHz) increment in the CMB intensity. The 
measured temperature decrement AFsz in the Raylegh- Jeans frequencies of the CMB is given by: 
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where T e is the temperature of the intra cluster medium (ICM), n e is the electron density, k% the Boltzmann constant, 
?cmb = 2.728°K is the temperature of the CMB, G T the Thompson cross section, m e the electron mass and f(v,T e ) 
accounts for frequency shift and relativistic corrections (c = 1 in our units). The integral is evaluated along the line of 
sight (l.o.s.). 

The hot gas in galaxy clusters is also responsible for X-ray emission due to bremsstrahlung and line radiation 
resulting from electron-ion collisions. The X-Ray surface brightness Sx is proportional to the integral along the line 




FIGURE 1. a) Contours in the q (1 — q\ plane for 38 galaxy clusters data [10] considering q(z) = qo + qiz- The best fit to the pair 
of free parameters is (qo,q\) = ( — 1.35,4.2). For comparison we have shown the straight lines denoting the transitions redshifts for 
two different flat ACDM models: z t = 0.9 for Q. A = 0.8 and z, = 0. 15 for Q. A = 0.43. b) Probability for the transition redshift. The 
best fit is Zi = 0.32. 



of sight of the square of the electron density: 
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where A e n is the X-Ray cooling function of the ICM in the cluster rest frame and Dl is the luminosity distance from 
the galaxy cluster. Observing the decrement of temperature of the CMB in the direction of galaxies clusters and also 
considering the X-rays observations, which are sensible to a different combination of the cluster electron density n e 
and the temperature T e , it is possible to break the degeneracy between concentration and temperature and to calculate 
the distance of angular diameter, yielding 

where 9 C refers to a characteristic scale of the cluster along the l.o.s., whose exact meaning depends on the density 
model adopted. This technique for measuring distances is completely independent of other techniques and it can be 
used to measure distances at high redshifts directly. 

Let us now consider the 38 measurements of angular diameter distances from galaxy clusters as obtained through 
SZE/X-ray method by Bonamente and coworkers [10]. The cluster plasma and dark matter distributions were analyzed 
using a hydrostatic equilibrium model (assuming spherical symmetry) that accounts for radial variations in density, 
temperature and abundance. 

In our analysis we use a maximum likelihood determined by a % 2 statistics 
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where Da j is the observational angular diameter distance, Od Ao , is the uncertainty in the individual distance, (7 sta t 
is the contribution of the statistical errors added in quadrature (w 20%) and the complete set of parameters is given 
by p = (Ho,qo,q\). The common errors are: SZE point sources ±8%, X-ray background±2%, Galactic N# < ±1%, 
±15% for cluster asphericity, ±8% kinetic SZ and for CMB anisotropy < ±2%. When we combine the errors in 
quadrature, we find that the typical error are of 20%. The systematic effects will not considered in this preliminary 
analysis. Systematic contributions are: SZ calibration ±8%, X-ray flux calibration ±5%, radio halos +3% and x-ray 
temperatute calibration ±7.5%. For the sake consistency, the Hubble parameter Hq has been fixed by its best fit value 
H* = mkm/s/Mpc. 




FIGURE 2. a) Contours in the qo — q\ plane for 38 galaxy clusters data [10] considering q(z) = qo + qiz/(l + z)- The best fit to 
the pair of free parameters is (qo,qi ) = (—1.43,6.18). b) Probability for the transition redshift and the associated best fit at Zt = 0.30. 
Comparing with Figs. 1(a) and 1(b) we see that the results are weekly dependent on the parameterizations. 

Linear Parameterization: q = qo + qiz 

In Figs. 1(a) and 1(b) we show, respectively, the contour in the plane qo —q\ (68.3%, 90% and 95.4% c.l.) and 
likelihood of the transition redshift from the Bonamente et al. sample. The confidence region (la) are —2.6 < qo < 
—0.25 and 13 < q\ < —3. Such results favor a Universe with recent acceleration (go < 0) and a previous decelerating 
stage (dq/dz > 0). From Fig. 1(a) we see that he best fits to the free parameters are qo = — 1 .35, = 4.2 while for 
the transition redshift is z t — 0.32 (see Fig. 1(b)). Note the presence of a forbidden region forming a trapezium. The 
horizontal line at the top is defined by q\ =0, which leads to an infinite (positive or negative) transition redshift. Note 
also that the segment at 45% defines the infinite future (z, = —1). In addition, one may conclude that the vertical 
segment on the left closing the trapezium is also unphysical since it is associated z, < — 1, thereby demonstrating that 
the hatched trapezium is actually a physically forbidden region (for a similar analysis involving luminosity distance 
see [15]). For comparison we have also indicated in Fig. 1(a) the transition redshifts Zt = 0.15 corresponding to a flat 
ACDM with £2 A = 0.43, as well as, z t — 0.9 corresponding to £1 A ~ 0.8. 



2nd Parameterization: q = qo + qn/l +z 

In Figures 2(a) and 2(b) we display the corresponding plots for the second parameterization. The confidence region 
(1 <t) is now defined by: —2.4 < qo < —0.5 and 13.5 < q\ < 0. Such results also favor a Universe with recent 
acceleration (qo < 0) and a previous decelerating stage (dq/dz > 0). As indicated in Fig. 2(a), the best fits to the 
free parameters are qo = — 1.43 and q\ =6.18 while for the transition redshift is a little smaller Zt = 0.3 (see Fig. 1(b)). 
It should be noticed the presence of the forbidden region (trapezium) with a minor difference in comparison with Fig. 
1(a), namely, as an effect of the parameterization, the horizontal line now is at the bottom. Note also that a decelerating 
Universe today (qo > 0) is only marginally compatible at 2(7 of statistical confidence. 

The results in the qo — q\ planes for both cases suggest that: (i) the Universe had an earlier decelerating stage 
(qi = dq/dz > for linear case and q\ > \qo\ for 2nd parameterization), and (ii) the Universe has been accelerating 
(qo < 0) since z ~ 0.3. A similar result has been previously obtained using SNe type la as standard candles by Shapiro 
and Turner[12]. 



PROSPECTS FOR PLANCK SATELLITE MISSION 

Let us now discuss the potentiality of the SZE/X-ray technique when future data from Planck satellite mission become 
available [19]. Planck satellite mission is a project from European Space Agency whose frequency channels were 
carefully chosen for measurements of thermal Sunyaev-Zeldovich effect. In principle, the Planck satellite will see 
(through SZE) about 30,000 galaxy clusters over the whole sky with significant fraction of clusters near or beyond 
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FIGURE 3. a) we have contours on the plane (q ,q\) from synthetic Planck data as defined in Table 1. In the pessimistic case, 
210 P (pessimistic projection) and 1050 O (optimistic projection) clusters were considered by assuming a random error between 
10% and 20% in the resulting angular diameter distances., b) we display the likelihood for the transition redshift in the pessimist 
projection. We obtain zt = 0.31 ± 0.03. 



redshift unity. However, since accurate angular diameter measurements require long SZE/X-ray integrations, we do 
not expect that all observed clusters might have useful distance measurements to constrain cosmological parameters. 
Therefore, it is interesting to simulate two realistic samples of angular diameter distances (ADD) by using a fiducial 
model to D'^ = D^z,-,^,^,//^ ), where H^, and q\ are the best fit values to the linear case obtained from 
Bonamente et al. sample[10]. 

Table 1 

Z range Clusters bins Clusters/bin ADDj Error 
(P, O) (P,0) (P,0) 

[0.0,0.5] 100,500 10 10,50 15%, 10% 
[0.5,1.0] 70,350 10 7,35 17%, 12% 
[1.0,1.5] 40,200 10 4,20 20%, 15% 

The first simulation (termed pessimistic - P), assumes that only 210 clusters are distributed in the redshift ranges in 
the following form: < z < 0.5 (100), 0.5 < z < 1 (70) and 1 < z < 1.5 (40) with ADD statistical errors of 15%, 17% 
and 20%, respectively (see Table 1). In the second one (optimistic case - O), 1050 clusters were redshift distributed 
as follows, < z < 0.5 (500), 0.5 < z < 1 (350) and 1 < z < 1.5 (200) with ADD statistical errors of 10%, 12% and 
15%, respectively. The redshift intervals were partitioned into bins (Az = 0.05) with the clusters distributed as shown 
in Table 1[20]. Both simulations were carried out by marginalizing over the Hq parameter in Da(z ( ,p) in Eqs. (1) and 
(5). 

In Fig. 3 we display the results of our simulations for the linear parameterization. The contours correspond 68%, 
95% and 99.7% c.l. for the optimistic (colored inner contours) and pessimistic case (outer contours), respectively. 
Comparing with Fig. 1(a) we see that the allowed region was remarkably reduced even in the pessimistic case. This 
means that ADD from SZE/X-ray will become a potent tool for measuring cosmological parameters fully independent 
and competitive with SNe type la. 

CONCLUSIONS 

We have shown that the combination of Sunyaev-Zeldovich/X-ray data from galaxy clusters is an interesting technique 
for accessing the present accelerating stage of the Universe. This result follows from a new kinematic approach based 
on the angular diameter distance of galaxy clusters obtained from SZE/X-ray measurements. By using two different 



parameterizations, it was found the existence of a transition from a decelerating to an accelerating expansion was 
relatively recent (z,t — 0.3). 

The ability of the future Planck satellite mission data to constrain the accelerating stage was discussed by simulating 
two realistic samples of angular diameters from clusters. The allowed regions in space parameter was significantly 
constrained for both the pessimistic and optimist simulations (Figure 3). The limits on the transition redshift derived 
here reinforces the extreme interest on the observational search for the combination of Sunyaev-Zeldovich and X-ray 
surface brightness from galaxy clusters. 

Finally, it should also be stressed that the present results depends neither on the validity of general relativity nor 
the matter-energy contents of the Universe and, perhaps, more important, it is also independent from SNe type la 
observations. 
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